%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Dielectric function for Aluminium.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;		addpath C:\Users\Philippe\Physik\Appendix\Appendix;		Constant;

Eps_infty=      1;
w_p      = 15.702;											%plasmon frequency
Gam_d    =   .696;
w        =[0:.1:40];										%incident light frequency

Eps1=1-w_p^2./(w.^2+Gam_d^2);
Eps2=w_p^2*Gam_d./w./(w.^2+Gam_d^2);

figure(1);				clf;
subplot(221),plot(w,0,'-k',w_p*[1 1],[-10 10],'-.k',w,1,'-k',w,Eps1,'r',w,Eps2,'g','LineWidth',1.4)
set(gca,'yLim',[-5 5],'FontSize',16,'LineWidth',1.4)%*E0/(2*pi*hP)
xlabel('\fontsize{16}{\it\omega}   (eV)');
ylabel('\fontsize{16}{\it\epsilon_r}   ,   {\it\epsilon_i}');
